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Abstract 

Casimir forces are conventionally computed by analyzing the effects of boundary conditions on 
a fluctuating quantum field. Although this analysis provides a clean and calculationally tractable 
idealization, it does not always accurately capture the characteristics of real materials, which cannot 
constrain the modes of the fluctuating field at all energies. We study the vacuum polarization energy 
of renormalizable, continuum quantum field theory in the presence of a background field, designed 
to impose a Dirichlet boundary condition in a particular limit. We show that in two and three 
space dimensions, as a background field becomes concentrated on the surface on which the Dirichlet 
boundary condition would eventually hold, the Casimir energy diverges. This result implies that 
the energy depends in detail on the properties of the material, which are not captured by the 
idealized boundary conditions. This divergence does not affect the force between rigid bodies, but 
it does invalidate calculations of Casimir stresses based on idealized boundary conditions. 
Keywords: Energy densities, Green's functions, renormalization, Casimir effect 
PACS: ll.lO.-z, ll.lO.Gh, 03.70. +k 
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1 Introduction 



The development of powerful experimental methods to measure the Casimir force and its 
potential importance in micromechanical devices have stimulated new interest in the phys- 
ical principles underlying the Casimir effect El El E] • The Casimir force was originally 
discovered by Casimir and Polder in their study of the long range, relativistic tail of the van 
der Waals force between neutral atoms [Sj. Immediately afterwards, Casimir showed that 
the force between grounded, conducting plates could be understood as a modification of 
the quantum fluctuations of the electromagnetic fields forced to obey conducting boundary 
conditions on the plates [H]. 

The intriguing idea of a force generated by the modification of the zero point energy of 
a fluctuating quantum field due to the imposition of a boundary condition has generated 
continuing theoretical interest in generalized "Casimir problems" [7j. The general idea is to 
specify a) an otherwise free quantum field (scalar, Dirac, gauge, etc.), b) a boundary con- 
dition (Dirichlet, confining, conducting, etc.), and c) a geometry (parallel plates, spherical, 
cyhndrical, etc.), and to calculate the forces, pressures, and other potential observables that 
arise when the field is forced to obey the boundary condition on the surface. 

These Casimir problems are idealizations in which the physical interactions between 
the fluctuating fields and matter have been replaced ab initio by boundary conditions. A 
real material cannot constrain modes of the field with wavelengths much smaller than the 
typical length scale of its interactions. The interactions become negligible at wavelengths 
less than certain physically determined cutoffs. In contrast, a boundary condition constrains 
all modes. To calculate the Casimir energy it is necessary to sum over the zero point energy 
of all modes. This sum is highly divergent in the ultraviolet and these divergences depend 
on the boundary conditions. Subtraction of the vacuum energy in the absence of boundaries 
removes only the worst divergence (quartic in three space dimensions). We ask whether the 
Casimir energy and other potential observables can be defined independent of the cutoffs 
that characterize the actual interactions between the fluctuating fields and the matter. If 
so, then one can define an abstract Casimir problem, depending on the field, the boundary 
condition, and the geometry alone. If not, then the Casimir energy is unavoidably entangled 
with the detailed material physics at hand. We study the simple case of a scalar field, (p, 
the Dirichlet condition, (p = 0, and three geometries: single plate, parallel plates and sphere. 
We believe this case illustrates general principles that affect all Casimir calculations. 

We find that the Casimir energy always depends on the cutoffs and diverges as the cutoffs 
go to infinity. The divergences in three dimensions are even more severe than the divergences 
in one and two dimensions described in earlier work jH]. We show that these divergences can 
be isolated in low-order Feynman diagrams, which can be studied with standard methods of 
renormalized quantum field theory. It is important to emphasize, however, that the Casimir 
force between rigid bodies, which is the only thing that has been measured experimentally. 
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is cutoff independent and therefore finite in the hmit of infinite cutoff. The most important 
example of an observable that is unavoidably cutoff dependent is the "Casimir pressure" 
{—dE/dA), which we study for the case of a sphere. 

The method of renormalization in continuum quantum field theory without boundaries 
(QFT) provides the only physical way to regulate, discuss, and eventually remove diver- 
gences. We propose to replace a boundary condition by a renormalizable coupling between 
the fluctuating field and a non-dynamical background field representing the material. When 
the background is smooth and the coupling is finite, then renormalization ensures that the 
zero point energy of the fluctuating field relative to the vacuum is finite. In the limits that 
the background becomes sharply peaked on a surface and the coupling becomes strong, all 
modes of the fluctuating field will obey a boundary condition on the surface. If a physical 
observable like the pressure diverges in this limit then the observable depends on the cut- 
offs that in reality keep the background from becoming arbitrarily sharp and the coupling 
from becoming arbitrarily strong. Our prescription is: a) fix the background, b) compute 
the renormalized zero point energy, c) take the "boundary condition limit". This is to be 
contrasted with the standard approach: a') take the boundary condition limit, b') compute 
the zero point energy. The two approaches disagree for all surfaces (plates and spheres) and 
in all dimensions that we have studied. We believe that our prescription for the order of 
limits is the physical one, and that the result obtained by imposing the boundary condition 
ab initio should be discarded when it disagrees with it. 

Here we study the imposition of Dirichlet boundary conditions on a scalar field. It is 
straightforward to write down a QFT describing a renormalizable interaction between the 
fluctuating scalar field and a static, non-dynamical background field 0"(x), 

Ant(0,(7) = -iA(7(x)0^(x,t) (1) 

and to choose a limit involving the shape of o"(x) and the coupling strength. A, between 
and a that produces the desired boundary conditions on specified surfaces. The interaction 
strength is characterized by two parameters: the width A of the background and the strength 
A. In the limit A — the background becomes a surface 5-function. We refer to this as the 
"sharp limit." Once the sharp limit has been taken, the limit A — > oo enforces the Dirichlet 
boundary condition on all modes of the fluctuating field. We refer to this as the "strong" 
limit. The boundary condition emerges on the surface only in the limit that the background 
becomes both sharp and strong. Both A and A represent physical cutoffs characteristic of 
the material with which the fluctuating field interacts. We identify A with the physical 
thickness of the surface and A plays a role similar to the plasma frequency: modes with 
frequency much larger than the scale determined by A are not constrained at the boundary. 

Ideally, we seek a Casimir energy that reflects only the effects of the boundary conditions 
and not any other features of o"(x). Therefore we do not specify any action for a except 
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Figure 1: The Feynman graph series for the Casimir energy. 



for the standard counterterms induced by the (f)-a interaction. The coefficients of the coun- 
terterms are fixed by renormahzation conditions apphed to perturbative Green's functions. 
Having been fixed in perturbation theory, the counterterms are fixed once and for all and 
must serve to remove the divergences that arise for any physically sensible o"(x). More- 
over, the renormahzation conditions are independent of the particular choice of background 
ct(x), so it makes sense to compare results for different choices of o"(x), i.e. different geome- 
tries. In the simple theory we are considering, the only counterterms are £cti(o") = Cicr(x), 
which renormalizes the tadpole graph and, in three dimensions, £ct2(c) = C20"^(x), which 
renormalizes the self-energy. 

The vacuum polarization energy — the Casimir energy, E[a] — for the field in the 
background a is given by the sum of all one-loop Feynman diagrams plus the contributions 
of the counterterms, as shown in Fig. ^ 

The diagrams generate a representation for E[a] as a power series in A. It is easy to show 
that 

• for piecewise continuous a (characterized by a width A), the contribution to the renor- 
malized Casimir energy at each order in A is finite, and 

• the sum over the terms in the power series in A gives a finite, renormalized Casimir 
energy for any fixed A and A. 

The first result is a general consequence of the renormalizability of the 0-cr field theory. 
The second result can be demonstrated explicitly by rewriting the Casimir energy as a 
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sum/integral over scattering data (bound state energies and scattering phase shifts) in the 
fashion of Refs. OEl^B' where it is manifestly finite. 

The thickness A sets the scale for the Fourier components, (t(p) = J d^x cr(x) exp(zp ■ x) 
in D space dimensions, of the background field , which are integrated in the one-loop diagrams 
of Fig. H As A ^ the Fourier integrals in low-order diagrams diverge. This simple 
mechanism is responsible for the cutoff dependence of the Casimir energy. Because it occurs 
in low-order diagrams (order < 3 for three dimensions or less) we can examine the divergences 
with well-known and universally accepted methods of analysis^n]- This analysis forms the 
technical core of this paper. 

The simplest diagram, the tadpole, is momentum independent. In our renormalization 
scheme it is canceled exactly by the counterterm i leaving no contribution to the Casimir 
energy at order A. The two-point function has a simple representation in terms of 5"(p). In 
two dimensions, we have 

and in three dimensions, after renormalization we have 

Note the dependence of -E'3q [(t] on the (spacelike) renormalization scale /i. In the first part 
of Section 3 we show that these contributions to the Casimir energy diverge like In A and 
(lnA)/A in the sharp limit. There is no mechanism available to cancel these divergences. 
In two dimensions, if we use the counterterm £ct2(o') with finite coefficient C2 to make a 
finite renormalization of the two-point function, the divergence only gets worse, as discussed 
in Section 3.2. In three dimensions we have already renormalized the theory as the fi de- 
pendence of Ej^ [a] attests and yet the divergence in the sharp limit persists. We show that 
these cutoff dependences affect the Casimir pressure. 

The reader may still worry that some more subtle interpretation of renormalization will 
still succeed in removing the divergence, especially in i?3j3 [(T] which had a primitive (loop) 
divergence prior to renormalization even for fixed A. To address that concern we study the 
three-point function in three dimensions, for which there is no loop divergence. We show 
that this diagram also diverges (like In A) in the sharp limit — the situation is analogous to 
the two-point function in two dimensions. There is no mechanism to remove this divergence 
at any stage in the formulation of the theory, and it contributes to the Casimir pressure as 
well. Thus the Casimir pressure on a sphere is unavoidably cutoff dependent both in two 
and three dimensions. 

We have made many consistency checks on this calculation. In addition to the analytical 
analysis of the A ^ limit given in Sections 3 and 4, we have computed (in Section 5) 
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the diagrams in the sharp limit with a (5-function source but with a frequency cutoff on 
a(p). This is a simpler computation and it yields the same result. Finally we have used the 
"interface" formalism developed in Ref. ^2] (in the Appendix) to generalize from one to two 
to three dimensions, again confirming the result of the direct analytical analysis. 

The fact that the energy of a fluctuating field diverges when a boundary condition is 
imposed has been known for many years f3[ E] • However divergences are common in quan- 
tum field theory and renormalization was invented to remove them. Special calculational 
methods have been developed to eliminate the extra divergences that arise in the presence 
of boundaries, which give finite, cutoff-independent Casimir energies j7j. These methods do 
not address the crucial question of whether the discarded divergences reflect physically sig- 
nificant cutoff dependence or are merely unphysical artifacts of the calculational method. 
Several studies of these divergences have been made. Symanzik studied the divergences of a 
scalar field theory in the presence of a surface interaction that imposes the boundary con- 
dition = ^3]. He showed that it is possible to introduce a set of surface counterterms 
sufficient to cancel all the divergences that arise in this theory. There are many independent 
counterterms and their forms are complex. However, the counterterms available in a renor- 
malizable quantum field theory are few in number and simple in form, merely polynomials 
in the field and its derivatives. There is no freedom to introduce Symanzik's ad hoc surface 
counterterms if one views the boundary condition as arising from an underlying renormal- 
izable quantum field theory (as it must). Therefore his methods cannot be used to explore 
the physical implications of the divergences generated by boundary conditions. 

The work of Candelas^l] and Candelas and Deutsch^H] is similar to our both in spirit 
and conclusions. In Ref. ^H] they studied the Casimir energy and energy density for scalar 
and electromagnetic fields in the presence of boundaries. Candelas later refined this work 
and extended it to boundaries between dielectrics. They point out that the vacuum polar- 
ization energy density generally diverges as one approaches a boundary and state that these 
divergences are not the ones that can be cancelled by standard renormalization methods. 
They note that the famous special cases where the divergences cancel - the conformal scalar 
field near a planar boundarv|16j and the electromagnetic field near a plane or a sphere |T7j - 
are indeed exceptional. The cancellations fail and the Casimir energy density is divergent if 
the surface has (arbitrarily small) shape imperfections. Candelas and Deutsch speculate that 
divergences are generic and that a more realistic treatment of the boundary would confirm 
this. In Ref. ^3] Candelas develops a model of a dielectric in order to examine the diver- 
gences in Casimir energies under more realistic conditions. He finds that divergences are 
generic. The divergences that arise in the dielectric case are worse than those found when 
conducting boundary conditions are applied ab initio. Candelas warns that his results cast 
doubt on Casimir calculations which begin with idealized boundary conditions. However, in 
the years following Candelas's work this warning seems largely to have been ignored. The 
dielectric model used by Candelas does not correspond to a renormalizable continuum quan- 
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turn field theory. ^ Therefore one cannot state for certain which divergences can or cannot 
be cancelled by counterterms. Also, when he studies conductors or the Dirichlet problem, 
Candelas imposes boundary conditions ab initio, introducing the possibility of missing phys- 
ically significant divergences. Our work can be regarded as a reformulation of Candelas's 
critique of Casimir calculations in the context of renormalizable quantum field theory. In 
this context there is no longer any reason to doubt the physical significance of the cutoff 
dependence we have found. 

Finally our results disagree with the recent results of Ref . fH] , which also considers low- 
order Feynman diagrams, but obtains a finite contribution to the stress in the boundary 
condition limit. The specific disagreements are discussed in the Conclusions after we have 
introduced the conceptual framework necessary for the discussion. 

We note that some of these results have appeared in a brief introduction to our work 
jS] and in a conference presentation where the one- dimensional case was discussed 
in complete detail. A summary of the two dimensional case was also given in [Hj. The 
divergences of the Casimir energy in one dimension are less severe, but still prevent one from 
defining a Casimir energy in the Dirichlet boundary condition limit. Finally, we note that 
subtleties similar to the ones we have considered were addressed in the context of dispersive 
media in Ref. [20] 



2 General structure of the Casimir energy 

In this section we outline the calculation of the Casimir energy of a scalar field in a piecewise 
continuous scalar background field in two or three space dimensions {D = 2,3). We review 
the Feynman diagram expansion and renormalization issues and summarize the scattering 
data approach which sums all Feynman diagrams and gives an unambiguous renormalized 
vacuum polarization energy. We take a first look at how divergences can arise in the renormal- 
ized Casimir energy as the background becomes singular and identify the Feynman diagrams 
that play a central role. 

We consider the theory of a real scalar field (j) coupled to a time-independent background 
field cr(x), described by the Lagrangian, 

C = -9,09^0 - - ^0V(x) + CcrW] . (4) 

£ct[c"] is the counterterm Lagrangian required by renormalization. 

The Casimir energy is defined as the vacuum energy in the presence of a minus the 
vacuum energy in its absence. It can be written formally as the sum over the shifts in 

^Candelas models the dielectric in terms of elementary electric and magnetic oscillators, which could in 
principle be embedded in a renormalizable theory. However the effective field theory obtained by integrating 
out the oscillators to obtain a diclcctic is no longer renormalizable. 
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the zero-point energies of all the modes of relative to the trivial background a = 0, 
-E'bare[c"] = ^ Sn('^n.[o"] — uJn^). Equivalently, using the effective action formalism, i^barcio"] 
is given by the sum of all one-loop Feynman diagrams with at least one external a field, as 
shown in Fig. Q The low-order diagrams generate divergences, which are canceled by the 
counterterms in Cqt- Combining its contribution to the energy, i?cT[o"], with ii^bare[c"] yields 
the renormalized Casimir energy E[a]. 

We have introduced cr(x) as a non-dynamical, external field, solely to constrain 0. In 
addition to the 0-cr coupling we include only the counterterms induced by the renormalization 
process. Of course one could always add non-trivial dynamics for a in the form C[a] = 
^{daY — V (a) . Such terms are undoubtedly present in a real material, where they correspond 
to the self-interactions of the material. However, since we want to know whether the vacuum 
energy associated with the fluctuations of can be isolated from the rest of the problem, we 
suppress any additional dynamics associated with a. 

Because we are calculating in a renormalizable quantum field theory, E[a] will be finite for 
any smooth a and finite A. Renormalization is straightforward. First the theory is regulated. 
Next, superficially divergent Feynman diagrams are identified. Only the one-point function 
(the "tadpole" ) and the two-point function (the a self-energy) are superficially divergent for 
D < 3. The divergences are local and can therefore be canceled by counterterms of the form 



Here is a regulator, for example the deviation from the physical spacetime dimension in 
dimensional regularization. The constants ci and C2 are fixed by renormalization conditions. 
We choose the "no tadpole" condition (a) = to fix ci, and we require that the renormalized 
two-point function, vanish at = —fi^ to fix C2- Of course any other choice of renormal- 
ization scheme can be related to ours by appropriate renormalization group transformations. 
We wish to stress that the renormalization process fixes the counterterms completely, inde- 
pendent of cr. So, for example, the counterterm ci is given by. 



in n spacetime dimensions. No ad hoc modifications of the counterterms can be introduced 
to cancel an unanticipated divergence. 

Although the Casimir energy can be formally represented as the sum of all one-loop 
Feynman diagrams, only the first few can be computed directly. As the number of external 
lines increases so does the number of integrations over the external momenta carried by 
5"(p), the Fourier transform of cr(x). In practice it becomes tedious to compute beyond the 
three-point function directly. Furthermore the power series in A generated by the Feynman 
diagrams may not converge. We rely instead on methods that allows us to compute the 



CctW] = ci{r])cr + C2(?7)(T^ 



(5) 




(6) 
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Casimir energy of the background configuration exactly for any A, while still performing all 
the necessary renormalization in the perturbative sector. The full Casimir energy can be 
computed exactly using phase shift methods, [HIEIE]) by summing the derivative expansion 
|211 1^ . or through analysis of Green's functions |1H I23j. All these methods show that once 
the counterterms corresponding to low-order diagrams have been included, the total Casimir 
energy E[a] for a smooth background is finite. Expanding this result in powers of A, we 
recover the diagrammatic expansion and can isolate the effects of the counterterms. 

Before analyzing the contributions of order A^ and A^ in detail, we give an overview of 
the possible divergences in E[a] when the background a becomes sharp and strong. We will 
consider a square barrier of width A and height 1/A. Any smooth shape that approaches a 
surface delta function as A — would do. 

• The tadpole diagram is local, i.e. proportional to cr(x), and therefore is canceled 
entirely by the counterterm Ci = Cicr(x). The resulting rule for calculation is merely 
to drop the contribution of first order in A to the Casimir energy. No matter how 
singular the source, the contribution of order A to vanishes as a result of our choice 
of renormalization scheme. 

• The contribution of order A^ comes entirely from the renormalized two-point function, 
which can be computed using standard methods. So the rule for calculation is to 
remove the terms of order A^ from the scattering data and replace them by -E'-^-'fcx]. 
The expressions for i?'^^)[cr] in two and three space dimensions have been given in 
eqs. (j21) and eq. Q respectively. 

^(2) 

[a] can diverge if the background a has a Fourier spectrum that vanishes too slowly 
at large p. This is precisely what happens when A ^ 0. Thus, examination of the 
renormalized two-point function alone is enough to establish that the renormalized 
Casimir energy diverges in the Dirichlet limit. 

• The higher order (in A) contributions to E[a] are guaranteed to be finite for smooth 
enough a and finite A because only the tadpole and self-energy diagrams have primitive 
divergences in the underlying field theory for D < 3. We find that cr(x) which are 
piecewise continuous are smooth enough to give a finite renormalized Casimir energy. 
Thus the renormalized Casimir energy is finite for any fixed A and A. 

In two dimensions there are no divergences beyond the two-point function even in the 
sharp limit. However, in three dimensions as the background field becomes sharp, i.e. 
as A — i> 0, the three-point function diverges as well. In Section 4 we calculate the 
three-point function in the background of parallel plates of width A and also in the 
background of a spherical shell of with A. We study the limit A — and show that 
this contribution to E diverges. 
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In Section 5, we calculate the contribution to the vacuum polarization energy at n 
order in A for a surface (5-function in three dimensions for both plane and spherical 
geometries. We find a result that is finite for n > 3, but diverges logarithmically for 
n = 3, a divergence that cannot have anything to do with choice of renormalization 
scheme or other subtleties. 

To summarize: Potential divergences in the Casimir energy in the Dirichlet limit are 
isolated in the two-point function in two dimensions and in the two and three-point functions 
in three dimensions. These contributions are analyzed in depth in the next two sections for 
single plate, parallel plate and spherical geometries in two and three dimensions. 



3 The Two-Point Function 



In this section and the next we analyze the low-order Feynman diagrams that are potentially 
divergent in the sharp limit, A ^ 0. We treat the case of a square barrier, although any 
one-parameter family of functions which approaches a (5-function as A ^ would do. In this 
section we study the renormalized two-point function as defined in eqs. (0) and Q. While 
this analysis alone is sufficient to establish a divergence in the total energy, the fact that 
in three space dimensions it requires renormalization may make some readers uneasy. The 
three-point function, while more difficult to compute, has no such complication and we will 
study its divergences in the next section. 
We study three geometries: a single plate 

a^{z) = ^{9iz + A/2) -9iz- A/2)), (7) 

parallel plates separated by a distance 2L (both configurations in D(< 3) space dimensions), 

a\\iz) = ^{e{\z\ -L + A/2) -ei\z\-L- A/2)) , (8) 
and a sphere of radius R in D=3 space dimensions, 

4. + 1)3 -g) (''('--«)-''(''-«-^))- P) 

We have chosen to normalize a analogously for plates and spheres. For plates we have 
normalized a so its integral across each plate is proportional to the surface area and inde- 
pendent of the width A and the number of dimensions, i.e. we require / dza{z) = N, where 
iV = 1 or 2, is the number of plates present. For a\ and ay we study S, the Casimir energy 
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per unit area in three dimensions or per unit length in two dimensions. In the case of a 
spherical shell of radius R we have normalized a so that 



/ 



(i^X(Jo(x) = 1 



(10) 



independent of the width A. We are aware of the fact that different normalization choices 
imply different dimensions for A. These can be easily introduced, for example, by introducing 
appropriate powers of m, the mass of the fluctuating particle. 

The choice of normalization has little significance for single or parallel plates. However 
it does have physical consequences when we measure the stress on a circular or spherical 
shell. To measure the stress, it is necessary to compare the energies of shells of radius R and 
R + SR. To make this comparison, one must decide how the matter behaves as R changes. 
Unfortunately, one of the shortcomings of the traditional boundary condition approach is 
that it tells us nothing about the physical properties of the material. So we must consider 
a range of possible behaviors. With the normalization chosen above, / drao{r) ~ 1/i?^ in 
three dimensions. This normalization corresponds to keeping the volume of a fixed as the 
shell is expanded, as if it were a physical substance. Another possibility would be to hold 
J drao{r) fixed as R varies. We consider both possibilities when discussing the stress on 
a shell. We will see that it is not possible to choose any normalization of a that would 
eliminate divergences in the renormalized Casimir stress on a sphere in three dimensions. 

It is convenient to work in momentum space when computing n-point functions. The 
background fields in momentum space are given by the Fourier transforms of eqs. (jH}, 
in the case of the single plate, parallel plates and sphere respectively: 



where q = |q|. 

Most of our interest lies in three dimensions where there has been some controversy about 
the nature of the divergences in this problem JHl- At the end of this section we present 
analogous results in two dimensions for comparison. The two-point function contribution to 
the Casimir energy can be computed analytically for a massless scalar field. Therefore we 
begin with the massless case. 





(12) 



3 qR [cos{qR) - (1 + A/R)cos{q{R + A))] - sm{qR) + sm{q{R + A)) 
^ (i? + A)3 -i?3 



(13) 
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3.1 Massless Case 



In two dimensions the two-point function contribution to the Casimir energy is infrared 
divergent in the massless case. So in this subsection we confine ourselves to three dimensions. 
In three dimensions the renormalized two-point function was given in eq. Q, and simplifies 
in the massless case to 

For the single plate, eq. ((Tj), we obtain 

^ 'ln(/iA)+7-l) , (15) 



327r2A 

for the energy per unit area, where 7 = 0.577 ... is Euler's constant. From this we conclude 
that the Casimir energy per unit area of a single plate diverges as ;^ In A in the sharp limit, 
A — s> 0. This is the first example of a general result: if we normalize the background field 
a to the surface area of the body under consideration, then the form of the leading and the 
first subleading divergence is always 

^S = -^^(ln(/iA)+7-l), (16) 

where A is the area. In the cases of single and parallel plates there are no further divergences 
in the sharp limit, but in the case of a sphere there is a further sub-sub leading divergence 
that is not universal. 

Next consider the case of parallel plates, eq. (jH)), 

^3dV||) = ^ / ^ - (1 + COS (2pL)) sin^ (pA/2) In (///x^) (17) 

This integral can be carried out by deforming the integration contour in the complex plane, 

= { In (A/i) + 7 - 1 + Q + ^) In (1 + A/2L) 



1 L 

2 ~ A 



ln(l - A/2L) 



It is easy to check that this result has the proper behavior as the plates separate to infinity or 
coalesce. As L ^ 00, the energy per unit area approaches twice the result for a single plate of 
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unit strength (see eq. (US))). And as the plates coalesce at L = A/2, S^^{a\\) becomes equal 
to the energy per unit area of a single plate of strength two with width 2A, as expected. As 
expected, the divergence in S^^{a\\) is twice that in S^^{a\). 

Although the Casimir energy for two parallel plates diverges as A ^ 0, the divergent 
terms are independent of L, explicitly, 

and therefore the two-point contribution to the Casimir force remains finite in the sharp 
limit. Also, only the L- independent terms depend on the renormalization scale /x. That is, 
this force does not have renormalization ambiguities. 

Next, we turn to the spherical shell. Using eq. Q, we obtain. 



p(2). , A2 (2 4(l + e)_ 2 + 26 

3° ^"""^ " 2567r3i?A2(l + e + '.e^y I 3 (2 + e)^ ''''' 2 + e 



-2e(ln(/iA)+7-l) + e2( -l + 21n^^ + 2(ln(/iA)+7-i; 



^n^±^ + ^(ln(/iA)+7-l) ) }■ (20) 



9 3 2 + e 3 

where e = A/R. The terms that diverge in the limit A ^ can be singled out 

E(g(ao) = "Y^^ ( 1^ (/^^) + 7 - 1) - ''^^) + fi^it^ for A ^ .(21) 

If we correct for the different normalizations, we find that the leading and the first subleading 
divergent contribution is the same (per unit area) as for the plate(s). However its physical 
consequences are quite different. In the case of an isolated sphere the stress, 

v'iM = -'-^ (22) 

is the quantity of interest. Apparently the stress diverges as -^l^lnA as A ^ 0. Had we 
normalized cTq in the same way as a\ and a\\, so that J drao{r) ~ constant, the R dependence 
of the prefactor in eq. (j2ip would have changed from oc 1/i?^ to oc i?^ . The stress changes 
sign, however it would still diverge as A ^ 0. This highlights the fact that the stress depends 
in detail on the surface dynamics. 

The reader might be tempted to eliminate the R dependence from the divergences in eq. 
fl21|) by choosing a different normalization for the background a and redefining R R+A/2. 
Such a transformation is possible, and would render the Casimir pressure from the two-point 
function finite. However, as we shall see below, the three-point function also contributes a 
divergent Casimir pressure, and no transformations can eliminate all of the divergences. 
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3.2 Case of Non-zero Mass 



Three Dimensions 

It has been claimed that the divergences in the sharp limit are an artifact of the specialization 
to massless bosons ^Hl- It is certainly true that low-order Feynman diagrams have infrared 
divergences in low dimension. However these are not the divergences that plague the Casimir 
energy. To show that the divergences of the previous section persist when m 7^ we present 
the case of a single plate. The other geometries of interest behave analogously. 
For a scalar boson of mass m in three dimensions we have to evaluate 



c(2)/ N ^' r' dp 1 . 2 
^3Dh) = T7-^/ ^^^^^ 




327r2A 




In (/iA) - 1 - — ^ / d^^^^ (1 -e-^"*^) - In(mA) 

(23) 




Note that the massless (m 0) and sharp (A 0) limits coincide in the ^ integral term 
and that , 

^^^3~^ (1 - e-«^) + Inr/j = -7. (24) 

In addition, the last term of eq. ()23p. the Feynman parameter integral, vanishes in the 

(2) 

massless limit. So the leading A ^ divergence of (c"|) is the same whether or not 
has a mass. This is true for all three geometries we have studied. Note that the subleading 
divergence proportional to is altered in the case of a massive scalar boson. 

When the mass of (p is non-zero it is possible to study the two dimensional case. This 
problem has been studied with a GauBian background in . For comparison with the three- 
dimensional case, we compute the Casimir energy of a single plate (line) in two dimensions 
with the background of eq. 

Two Dimensions 

In two spatial dimensions the vacuum polarization diagram is finite and requires no renor- 
malization. However, the contribution to the energy becomes infinite in the sharp limit, as 
we will now show. The integration in eq. Q cannot be performed analytically for arbitrary 
A. However the leading behavior as A — > can be easily obtained. The energy per unit 
length is given by 

S^^{(Ti) = In (Am) + terms finite for A 0. (25) 

327r 
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This calculation demonstrates that the contribution to the energy from the second-order 
Feynman diagram in two spatial dimensions is divergent in the limit A — 0, and this 
divergence cannot be renormalized away. If we try to renormalize the vacuum polarization 
graph (by a finite counterterm proportional to a^), the divergence only gets worse. The 
two-point function 

n(|p|) = -^tan-iM (26) 
47r|p| 2m 

is renormalized to 

n,en(|p|) = n(|p|)-n(/i) (27) 

where /i is a renormalization scale and the leading divergence now is proportional to 1/A. 
A numerical analysis of the vacuum polarization energy for a GauBian background peaked 
on a circle of radius R can be found in Ref. jTTj. There the energy density is studied, and 
the log-divergence that develops as A ^ can be seen to be concentrated on the emerging 
surface at R. 



3.3 Comments 

At this point already some important conclusions may be drawn: first, the renormalized 
second order Feynman diagram diverges in the sharp limit for all the simple geometries we 
have considered. Second, since the Feynman diagram expansion is an expansion in powers of 
A, the contribution from the second order diagram cannot be canceled by the other diagrams 
in the series. They are of different powers in A. Third, the diagram has been renormalized, 
and the counterterm has been fixed. Thus, the divergence in the sharp limit is physical 
and continuum quantum field theory does not provide a way to get rid of it. If one is only 
interested in forces that do not require the body to be deformed, one need not be troubled 
by cutoff dependence in the energy. However, if one has to vary the surface size or geometry 
as in the case of a stress, the sharp limit gives a divergent result. It is already clear at this 
point that throwing away divergent terms on an ad hoc basis will eliminate parts of the 
physical measurable stress; therefore, in cases where one has to deform the body to obtain 
the stress the sharp limit cannot be taken. 



4 The Three-Point Function 



To further clarify the nature of the divergences in the Dirichlet Casimir energy, we consider 
the three-point function in a scalar theory in three dimensions. No divergences originate 
from the loop integrations in the Feynman diagrams at this order (and beyond). As a result. 
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we do not need to discuss the subtle issue of regularization and renormalization and no 
divergence in the three-point function can be renormahzed away. It must instead reflect 
a physical cutoff dependence generated by the interaction of the fluctuating field with the 
background matter. 

The three-point function for a fluctuating boson field reduces to 

-(3)r,i _ ^-^^ f ^^^2 d^qs cri{qi-q2) 5-4(^2 - gs) 5-4(^3-^1) 



3D ^^i - ^ J (27r)4 (27r)4 (27r)4 - m'^ + ir] ql - m'^ + irj qf -m? + i7]' ^ ' 

Here a^iq) = J d'^xa{x)e^'^'^ = 27r6{q^)a{q) is the Fourier transform of the time independent 
background field (t{x) in four dimensions. The factor of ^, where T is a large time interval, 
converts the effective action for arbitrary a to the Casimir energy for time independent a. 

We consider two cases for o"(x): First we assume a spherically symmetric barrier of radius 
R and width A, given by eq. (0), and second, we consider parallel plates of width A separated 
by 2L, given by eq. (jHl). 

As already noted, the integral over the loop momentum in eq. ()28p converges. Thus 



-^3D W] evaluated without reference to regularization. It can be cast into the form 

.(3)r 



E',-^[a]^ / dPfiP) (29) 
Jo 

where it is understood that -EsDicr] is to be replaced by the energy per unit area, i^3D['^]) 
in the case of parallel plates. In eq. ()29|) and below we omit constants of proportionality 
that are independent of the parameters {R, L, A) characterizing the background field, which 
are irrelevant for the present discussion. For simplicity of notation we also suppress the 
dependence of E and / on R, L, A, A, etc. 

For the spherically symmetric background we find the integrand 

(P) = p4 r ^0 gij^2 0^^(^p 0^ 0-^ 



j dz (To{PV1 + z sin 26) g{P cos 9, P sin 6, z) 



where 



dx 



^jAiv? + 4x(l — x)q'^ + {p + 2xqzy 
X < tanh 



_^ 2xp — p + 2xqz 



Am^ + 4x(l — x)g^ + (p + 2xqzy 

. tanh- ^ + '"^" \ . (30) 

y 4m2 -|- 4x(l — x)q'^ + {p + 2xqzY 
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■ 1 2 3 4 5 6 7 

ln(P) 

Figure 2: The logarithm of the integrand, eq. of the three-point function for the spherical 
background, eq. ifl^) for various values of the thickness, e = A/R, of the shell. Energies are 
measured in m, the mass of the particle in the loop. For orientation we have also included a 
limiting function, f{P) ~ 1/P, for which E^^^ diverges. Here the case R = 0.8/m is shown. 

The integrals in eqs. ()29|1 and (p?IH) cannot be evaluated analytically so we study them nu- 
merically. 

In Fig. 121 we display a typical example for the integrand of the spherical problem, eq ()30|) . 
For a given non-zero width, A, there is always an intermediate region in which fo{P) oc 1/P. 
Beyond that region the integrand drops off faster than 1/P so the integral converges for 
finite width A. As the width decreases, the region in which /o(P) oc 1/P increases and 
continues to larger values in P. As A ^ 0, this behavior persists for arbitrarily large P and 
therefore the integral eq (j29p diverges in that limit. The upper boundary of the region that 
exhibits the 1/P behavior increases approximately like 1/A. We can therefore estimate the 
divergence as a function of the width, 

poo nPo /-l/A 

/ dPUP)^ dPfo{P)+ dP- = C\nA + ... (31) 

Jo Jo JPo 

where C is a constant and Pq denotes the momentum at which the 1/P behavior sets in. 
This logarithmic divergence is expected from the discussion above: in essence, 1/A serves 
as a large momentum cutoff on /o(P)- 

This divergence is studied in Fig. IHl where we show E^jy{ao) as function of the width 
parameter for various values of the position of the spherical barrier. Over the considered 
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2.50 




1.00 1 ' ^ ' ^ ' ^ ' 1 

0.00 0.05 0.10 0.15 0.20 

e 

Figure 3: The contribution of the three-point function to the energy as a function ofe = A/R. 
Energies are measured in units of the mass m of the particle in the loop. In those units the 
coupling constant X is chosen to be 1 and a scale is chosen to stress the R-dependence as 
well as the logarithmic dependence on A. 



range, A varies by about two orders of magnitude. On that scale the curves shown Fig. |21 
are approximately constant and thus verify the logarithmic divergence estimated in eq. ()31|) . 
So the divergence is not as severe as for the renormalized two-point function, but it is still 
present. As illustrated in the next section, where we discuss the (5-function background 
{ie. A ^ 0), this is the expected behavior. When we restore the dimensional factors we 
see that E^{ao) oc -^InA, so the constant of proportionality is R dependent. Therefore 

the stress ( —dE^^{ao)/dR ) also diverges logarithmically as A 0. As discussed in the 
previous section, the exact form of the R dependence depends on the normalization we 
choose for cxo. If we switch to the normalization used for plates where J dra{r) = 1, then 
£■30 (cTo) oc R^ In A, so the divergent contribution to the stress, and even its sign(!), depends 
on the physical interpretation we place on a, which is by no means clear. Note that the R 
dependence of this divergence differs from that of the two-point function, so for every choice 
of normalization of a the sum of the two terms diverges. 

The case in which the barriers approximate two parallel plates is simpler because one 
more integral can be evaluated analytically. The corresponding integrand in eq. ()29|) then 
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reads, 



/ll {P)=P I de S(P cos e, P sin 9) g{P cos 9, P sin ( 
Jo 



where 



q) = (T|j (p)(T|| {q)a\\ {p + q) , 
= ^ cos(gL) sin 



and 



dx 



pJo ^4m2 + 4x(l - x)g2 + (p + 2xg)2 

f , _i 2a;p — p + 2xq 

X < tanh 



Ay4m2 + Ax{l — + (p + 2xq)'^ 
-tanh-^ ^ ^ + '^^ |. (32) 

In Table^we present numerical results for the energy per unit area, £^3^ (cr||), which originates 
from the three-point function. Like the sphere, £^3^ diverges like InA as A ^ 0. However, 
this time the divergence is independent of L. As we found for the two-point function, the 

force between parallel plates due to the three-point function is finite and unambiguous as 
A ^ 0. 

To summarize, we have found that the contribution of the three-point function to the 
Casimir energy of a sphere or the Casimir energy per unit area of parallel plates diverges 
logarithmically in the sharp limit in three dimensions. We expect that this divergence, 
like that of the two-point function, is the same for all shapes. This is because the width 
parameter essentially plays the role of a cutoff on the Fourier spectrum of the background, 
independent of its geometry. There is no counterterm in the continuum field theory to 
cancel this divergence. The stress on a sphere suffers the same divergence. However the 
force between parallel plates does not. Instead it reaches a finite limit as A ^ 0. In the 
case of the sphere, the R dependence of cro(r) could be chosen to render the divergence 
independent of R. However no choice can remove the divergences in both the two and 
three-point functions. We conclude that the Dirichlet-Casimir stress on a sphere in three 
dimensions depends in an unavoidable way on the cutoffs that characterize the material. It 
is not possible to compute even the cutoff dependent term in the stress on a sphere, to say 
nothing of the finite parts, because it depends dramatically on what is assumed about the 
R dependence of o"o(r). 
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A 


L = 0.5 


L = 0.8 


L = 1.0 


L = 1.2 


L = 1.5 


0.010 


-163.04 


-160.24 


-160.88 


-161.04 


-161.28 


0.025 


-122.16 


-119.68 


-119.68 


-119.76 


-120.00 


0.050 


-95.36 


-92.64 


-92.16 


-92.08 


-92.24 


0.100 


-70.56 


-67.84 


-67.36 


-67.28 


-67.20 


0.200 


-48.08 


-44.72 


-45.12 


-45.04 


-44.96 



Table 1: The contribution of the three-point function to the energy per unit area, 
as a function of the width, A, and the distance of the parallel plates, L. Units are set by 
the mass, m, of the particle in the loop. Constants of proportionality, that do not vary 
with either A or L are omitted. The variation of £^^\a\\) across the columns of the Table is 
consistent within the accuracy of the calculation, with the divergent part of E^^^a^^) being 
independent of L. 

5 Vacuum Polarization Energy in a Sharp Background 



In this section we take a different approach to the calculation of the Casimir energy in a 
sharp background in three space dimensions. Instead of introducing a finite thickness, A, as 
a cutoff, we start with a 5-function background, a*, from the start and introduce a cutoff, 
A, on its Fourier spectrum. This enables us to compute all n-point functions for the parallel 
plate and spherical geometries and to explore the divergences in a simpler mathematical 
context. 

We first consider the case of parallel plates. To begin we impose the constraint that a* 
depends only on X3 and is symmetric under reflection in this coordinate, 

W = ^11(^3) , cr*\ixs) = a-||(-X3) . (33) 
The appropriate form of the free Green's function for this problem reads 

^ia.(x-i') r dp C0S{X3P + f) COs{x'sP + ^f) 

e=o,i 



J (27r)3 J_^ 271 f-f -p^ -m^ + ir] 



Here = {vo,Vi,V2) comprises the time component as well as those spatial components of 
a four vector, v^, orthogonal to X3. The quantum number £ = 0, 1 is associated with parity. 
The 77,-point function, T^'"'\a*^), is given by 

r(")(cr||) = J d'^Xi J d^X2...J (i^X„G(Xi,X2)cr||(x2)G'(x2, X3)cr||(x3) . . .(^(Xn,^!) Cr||(xi) 

(35) 
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We have already established that T^'^\a*i^) diverges for n = 1,2, and 3. They should be 
regarded as regulated by a cutoff which will be introduced explicitly below. For a background 
of the structure given in eq (jH^ the integrals / cPxi are straightforward (x = {x^ , , x"^)) . 
They yield overall (infinite) factors of the "time interval" and "area" of the plates. The 
energy per unit area is obtained by dividing out these factors. The order n contribution is 
(after Wick rotation) 



(-A)" f 



2n 



(2vr)= 



dpi 



27r + pi + m 



2 ■ ■ ■ 



dpr, 



27r q"^ + Pn + rn? 



dZi... dZnt{Pn] Zn, Zi)cr|| (Zi )t (p„; Zi, 2;2)cr|| (^2) . . . t{pn] Z^-l, Zn)(T\\{Zn) , 

(36) 



where 



t{p] z, z) = cos{zp + — ) cos{z'p + — ) 



2' ^ ' 2 

This expression simplifies dramatically when we assume a (5-function background 

al{z) = 6{z - L) + 6{z + L) (3. 
The integrals over Zi not only become not only trivial, but identical, thus 

dp cos{pL + 



(37) 



4SH,A) 



-2A)" r d^q 
2n J (27r)3 



1 

"2n 



d'^q 
(2^ 



E 

' €=0,1 

-A 



in-' 
2 ' 



27r + p2 + 



Ml) 



{(1 + 



-2uj{q)L\n 



r + (i 



-2Lu{q)L\n 



r} , (39) 



where the cutoff A appears as an argument for S to remind us that this is a regulated 
Casimir energy. We would like to mention that there is no physical interpretation of this 
cutoff because it cannot be directly linked to either the renormalization scale or the Fourier 
spectrum of the background field. This cutoff merely serves to render eq. (jH^ finite for all n. 
Simple power counting indicates that the n = 1, 2, and 3 contributions diverge as A — > cx), 
while, as expected, the contributions of order n = 4 and higher are finite. 
Keeping the cutoff fixed, all orders are finite and can be summed. 



d^q 



In 



(27r)3 



A 



+ 



A^ 



uj{q) Auj^iq) 



(1 



-4uj{q)L 



(27r) 



-t Vt"^ — m? In 



1 + I + 4J5(1 



-4^L^ 



(40) 
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It is easy to verify that the divergences of £^3D(o"y,A) as A ^ oo correspond to those we 
analyzed diagram by diagram earher. However the divergences are L- independent, so that 
we can extract a force per unit area (pressure) between the two plates, 



e 



^3dK,a) = -- — ^^ = -^ dt ^ , ; , — — -. (41) 



2 dL STT^y^ i + A + 



Q-itL\ 



which is finite as A ^ oo. Furthermore, it remains finite in the strong limit in which A ^ oo 
and the situation reduces to two Dirichlet plates. We then obtain the established result for 
the force between two Dirichlet plates as a function of m and In the Appendix we 

show that the same result is obtained in the framework of the interface formalism of Ref . ^2] • 
Next we consider the analogous calculation for a sharp background with spherical sym- 
metry, cr*(x) = 6{r — R)/4ttR'^. The starting point is the spherical decomposition of the free 
Green's function 

where ji are spherical Bessel functions and Y^m are spherical harmonics. We substitute this 
decomposition into the expression for the n-point function, eq. Again the integral over 
the time and angular coordinates is trivial. It can be factorized to extract the n-th order 
contribution to the energy E^'^{a*), 



k^dkn 



u)"^ + k'i + rn^ 
t - - - " 

POO 

drirlji{knri)a*{ri)je{kiri)... / drnrlji{kn-^irn)o-*{rn)je{knrn) (43) 



where the frequency integral has been been Wick-rotated. As in the case of parallel plates, 
the contributions for n = l,2,and3 diverge. They should be regarded as regulated by a 
cutoff to be introduced explicitly below. 

For the 5-function background the momentum integrals can be done, giving 

■ r roo j^(/ti?)^ 



- 1 ( ~^ 



uj'^ + + rri 



( 4^ ) Jo ^ [^'^'/' (i?vW^) [rVi:^^t^ 



(44) 



where and Ki, are modified Bessel functions. The sum over i is finite for fixed u, but 
divergences appear in the u integral for small n. To regulate these divergences we have 
introduced a cutoff. A, on the frequency integral. 
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-10.0 



U=2000 



1.5 2.0 2.5 1 . . 3.0 3.5 4.0 

ln(z) 

Figure 4: The integrand for in eq. diij) as a function of tlie integral variable z as given 
in eq. Tlie maximal angular momentum used in tliis numerical study is ima^ = 2000, 

which is sufficient in the uj range considered. 



The central question is the behavior of the frequency integral after summing over orbital 
angular momentum i. For n = 1,2 we expect divergences because we know that the loop 
integral in the Feynman graph diverges. Since we did not distinguish between loop and 
external momenta, Feynman loop integral divergences will show up as divergences in the 
frequency integral. For n = 1 we know that YliC^^ + ^)je{z)'^ = 1- After summing i, the 
k integral diverges linearly, as does the u integral. These two divergences combine to a 
quadratically divergent object as expected for the tadpole graph. When we renormalize, the 
entire n = 1 term is canceled by a counterterm. In figure lU we display the integrand, 

^max 

=5^(2^ + 1) [li+i/2{z)K,+y2{z)Y , z = RVu;^ + m\ (45) 

for n = 2,3,4 in a log-log plot. The sum over £ is performed numerically up to ^max large 
enough to ensure convergence at the values of u considered. This numerical study confirms 
the results we obtained by studying the n = 2 and n = 3 cases in Sections 3 and 4. For 
n = 2 the integrand goes like u;""^ with < a2 < 1. The best fit to a2 is 0.2064. So 
the divergences of the frequency integral is stronger than logarithmic, i. e. stronger than 
the loop divergence. For n = 3 our best fit to 0^3 is 1.0004, consistent with the logarithmic 
divergence that we expect from the analysis of Section 4. Finally, for n = A 1^ behaves like 
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l/cu^ (the best fit to is 1.9829) yielding a finite frequency integral. So do the higher orders 
in n. This confirms our earlier results that for a ^-function background the singularities are 
still contained in the lower order Feynman diagrams and are stronger than those from the 
loop integrals. 

Finally, we stress that unlike the case of parallel plates, the divergences in the case of 
a sphere do depend on the radius of the sphere and therefore afflict the Casimir pressure, 
making it cutoff dependent. 



6 Conclusions 

We have seen that the Casimir energy of a Dirichlet boundary can be modeled as the limit 
of the Casimir energy of a background potential, which can then be calculated using the 
conventional tools of renormalized quantum field theory. Although this energy is finite and 
unambiguous for any particular smooth background, the limit of this process can still diverge, 
implying that the energy depends in detail on the properties of the physical material used 
to impose the boundary conditions. For three space dimensions, the Casimir energy has 
In A divergences both for parallel plates and for the sphere. These divergences should not 
be confused with the loop divergences of Feynman diagrams: they also appear in diagrams 
that have no quantum field theory divergences. In particular, the third order diagram needs 
no renormalization but also diverges in the sharp limit, going like In A. 

This cutoff dependence may or may not enter into physically measurable quantities. The 
force between rigid bodies is never affected, so our analysis does not alter standard results 
such as the force between parallel plates. However, it does render meaningless calculations 
of stresses, such as the Casimir surface tension of the sphere. Such quantities are cutoff 
dependent, and cannot be defined independently of the material properties that determine 
the cutoff. 

Our results disagree with those of Ref. JHI; which also considers low-order Feynman 
diagrams, but obtains a finite contribution to the stress in the boundary condition limit. We 
disagree with this calculation on two grounds: 

• In Ref. PH] the background fields are analytically continued to fractional dimension. 
In dimensional regularization, one continues loop momenta to fractional dimensions, 
but background fields, which are simply complex-valued functions, remain fixed in 
the physical dimension (see for example j^HJ; P- 662). Analytically continuing the 
background fields introduces a dependence into the counterterm coefficients ci and C2 
(see eq. (0)). Counterterm insertions in other diagrams would then render the theory 
unrenormalizable . 
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• In Ref. ^H] it is claimed that the contribution of the two-point function to the Casimir 
energy is finite for D = 3 space dimensions because the divergences appear in the form 
r(l — -D/2), which has poles at D = 2 and D = 4 but not D = 3. However, this is a 
familiar misinterpretation of dimensional regularization. For example, divergences in 
QED in d spacetime dimensions appear as r(2 — d/2), but QED in 4 + 1 dimensions 
is not a finite theory. Similarly, the divergences of 0^ theory come from r(l — d/2), 
but this theory is not finite in 2 + 1 or 4 + 1 dimensions. The proper prescription in 
dimensional regularization is to use counterterms to cancel all divergences appearing 
in dimensions less than or equal to the physical dimension. 

If the results described here generalize to the electromagnetic case, as we expect they do, 
they would invalidate Boyer's result [T7j that the conducting sphere experiences a cutoff- 
independent, repulsive Casimir stress. For such a result to be correct, it is necessary to 
show that it can be obtained as the limit of an underlying smooth, renormalizable quantum 
field theory. Otherwise, the Boyer problem cannot be studied without reference to material 
properties, and the assertion that sphere has repulsive Casimir force is unwarranted. 
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7 Appendix: Interface Formalism and Force between 
Plates 

In this section we will verify that the interface formalism of Ref. ^2] applied to a (5-function 
type background gives the standard result for the force [25 between two parallel Dirichlet 
plates. 

The parallel plate geometry is a standard example of an interface problem, where n 
"trivial" spatial dimensions parallel to the plates are irrelevant for the QFT dynamics and 
the equations of motion become effectively m-dimensional. In the particular case of the 
parallel plates in three dimensions, m = 1 and n = 2. 

From the point of view of scattering theory, the interface is characterized by a separation 
of the wave functions ^(x) ~ ip{x±)e^^'^^^ where x = x± + x\\, leading to a continuous 
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degeneracy of bound and scattering states. The phase shifts 5{k) will only depend on the 
m-dimensional momenta k in the non-trivial directions, while the integration over the n- 
dimensional trivial momenta p leads to an apparent divergence in the total energy. In 
Ref. [12], it was shown that this spurious divergence is actually absent due to identities in 
scattering theory that generalize Levinson's theorem [23]. It is then possible to renormalize 
the Casimir energy per unit area of the interface, £ = i?//lpiato, with just the standard field 
theory counterterms. For the interesting case m = 1, n = 2, the formula for S reads [23 



127r ^ 



+ 



C-(2) , r- 
^FD + ^CT 



h{k) = {e+fi'')^ -^efi-fl\ (46) 

where /i is the mass of the fluctuating boson, /i^ — k| ^ = Uj^^ are the bound state energies and 

the counterterm for the first order diagram Sp^ has been chosen to cancel it completely. The 
specific form of h{k) originates from integrating out the n trivial momenta in dimensional 
regularization. 6e{k) is the phase shift in the i^^ partial wave, and 6^^^ and are its 
first and second Born approximations. Since we are working in three spatial dimensions, a 
second Born subtraction has to be performed in the phase shift 5^{k), and the corresponding 
second order diagram £^-^ has to be added back in. This diagram is then renormalized by 



a standard counterterm £ct- The final result, eq. ()46|) . yields a finite expression £ for any 
smooth background a{xj_). 

The phase shifts in eq. are defined through the m = 1 dimensional problem of a 
massive boson scattering off a potential cr{x) on a line, and £ = ± denotes the two parity 
channels of that problem. We denote the non-trivial coordinate perpendicular to the plates 
by X and write x = x + x\\ = x ej_ + Z/ieji"^'' + y2&i^'^ ■ 

Before we study the sharp limit a\\{x) = 6{x — L) + 6{x + L), it is useful to follow [12] 
and rewrite eq. (|46|) as an integral over the imaginary axis k = it using analytic properties 
of scattering data. This approach combines the scattering and bound state contribution 
to eq. ()46|) in a simple and concise formula which is also better from a numerical point of 
view 

The general procedure is explained thoroughly in Ref. ^21- Here, it suffices to mention 
that the phase shift may be expressed in terms of the Jost functions, 25^{k) = ^ In F^{k) / Fi{—k). 
Since the integrand in eq. is odd in /c, we can extend the integral to [— oo, oo] and de- 
form the contour to a large circle in the upper complex fc-plane, where Fi{k) is meromorphic. 
Owing to the two Born subtractions, the integrand falls fast enough at |fc| — >■ oo. In addition, 
the residues from simple poles at the bound states k = iKj cancel the explicit bound state 
contribution in eq. pUj) . leaving only the discontinuity across the branch cut of the power 
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(fc^ + m^) 2 in h{k), 

oo 

^ = ^ f dttVW^Yl h(^) - ^^^'^(^) - '^^'^(^)] + 42 + ^CT (47) 

where ^'^(t) = lnF^(A; = zt) is a real function of t. 

Now we are prepared to study the sharp hmit cr{x) = 6{x + L) + 6{x — L), which we 
interpreted as a two-channel problem on the half-line x > 0. Using the well-known Jost 
functions in the symmetric and antisymmetric channel, respectively, 



F^(k) = 1 + 1^ \l± e^'^^] =^ z/±(t) = In 
2k 



A(l±e-2^*) 



(4J 



it is readily seen that the t-integral in eq. (07j) diverges even at finite coupling A. This is the 
divergence of the Casimir energy in a sharp background from our previous discussion. The 
t-integral in eq. (j47|) is a compact way of summing all Feynman diagrams of order three and 
higher. As shown in Section 5 the divergence of this sum in the sharp limit is exclusively due 
to the third order diagram £^-^. To see this, we subtract one more Born term in eq. (P7j) and 
add back in the corresponding Feynman diagram Eyt)- Then the t-integral (now the sum of 
all diagrams of order four and higher) will be finite, while the third order diagram (which 
would be finite for a smooth background) diverges ^, c/. eq. (fH^ . 

oo 

4S-^/^*^^[l + 3e-]. (49) 



487r2 j t2 

m 

There is no field theoretic counterterm to renormalize this diagram and thus the divergence 
of the Casimir energy implies a real physical cutoff dependence. 

In the present case, however, this cutoff-dependence has no physical consequences since 
the logarithmic divergence in eq. (|49|) comes from the first term in the brackets, which 
turns out to be independent of the distance d = 2L between the plates. The same holds 
for the second order diagram Sp^). As a consequence, the physically measurable pressure 
V = —d£/d{2L) on the plates is finite even though the total Casimir energy diverges. This 
is again just our main statement about the cutoff independence of the Casimir force between 
rigid bodies. 

After a little algebra, the subsequent strong coupling limit A ^ oo yields a finite result^ 



^Note that it is (far) simpler to compute this diagram by (Born) expanding v^it) in powers of the couphng 
A and inserting in eq. (|47|l . than to compute the Feynman diagram directly. 

■^Note that the finite contributions to the force cancel between the Feynman diagram and the correspond- 
ing Born subtraction for both the 2'^'* and 3'^'^ order contributions. This is important since those terms, 
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for the Casimir pressure between two parallel plates in three dimensions. 



where = 2L is the distance between the plates. This result can also be derived from eq (j^T| 
with A ^ oo and A oo. The integral C{s) becomes particularly simple in the massless 
case, where we recover the standard result 

P(d) = -, m = 0. (51) 

^ ^ 480 ' ^ ' 
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